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Abstract 

In the usual Hofava-Lifshitz cosmological models, the scalar field is responsi- 
ble for dark matter. Using an additional scalar field, Saridakis [1] has formulated 
Hofava-Lifshitz cosmology with an effective dark energy sector. In the paper [1] the 
scalar fields do not interact with each other, here we extend this work to the inter- 
acting case, where matter scalar field <j) interact with dark energy scalar field a. We 
will show that in contrast with [1], where <7-filed is absent, we can obtain vf^ < — 1, 
that is we result to an effective dark energy presenting phantom behaviour. This 
behaviour is pure effect of the interaction. 
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1 Introduction 



Recent observations from type la supernovae [2] associated with Large Scale Structure [3] 
and Cosmic Microwave Background anisotropies [4] have provided main evidence for the 
cosmic acceleration. The combined analysis of cosmological observations suggests that 
the universe consists of about 70% dark energy, 30% dust matter (cold dark matter plus 
baryons), and negligible radiation. Although the nature and origin of dark energy are 
unknown, we still can propose some candidates to describe it, namely since we do not 
know where this dark energy comes from, and how to compute it from the first principles, 
we search for phenomenological models. The astronomical observations will then select 
one of these models. The most obvious theoretical candidate of dark energy is the cos- 
mological constant A (or vacuum energy) [5, 6] which has the equation of state parameter 
w — — 1. However, as it is well known, there are two difficulties that arise from the cosmo- 
logical constant scenario, namely the two famous cosmological constant problems — the 
"fine-tuning" problem and the "cosmic coincidence" problem [7] . An alternative proposal 
for dark energy is the dynamical dark energy scenario. This dynamical proposal is often 
realized by some scalar field mechanism which suggests that the specific energy form with 
negative pressure is provided by a scalar field evolving down a proper potential. So far, 
a plethora of scalar-field dark energy models have been studied, including quintessence 
[8], K-essence [9], tachyon [10], phantom [11] and quintom [12], and so forth. An alter- 
native way of explaining the observed acceleration of the late universe is to modify the 
gravitational theory and in the simplest case replace R with f(R) in the action which 
is well known as f(R) gravity. Here f(R) is an arbitrary function of scalar curvature 
(for recent reviews see [13, 14]). However, most of /(i?) -gravity models do not manage 
to pass the observational and theoretical tests (solar system, neutron stars and binary 
pulsar constraints), giving also rise to an unusual matter dominated epoch and leading 
to significant fine-tunings [15]. 

Recently Horava proposed a renormalizable gravity theory with higher spatial derivatives 
in four dimensions which reduces to Einstein gravity with a non-vanishing cosmological 
constant in IR but with improved UV behaviors [16, 17, 18]. It is similar to a scalar field 
theory of Lifshitz [19] in which the time dimension has weight 3 if a space dimension has 
weight 1, thus this theory is called Horava-Lifshitz gravity. This theory is not invariant 
under the full diffeomorphism group of general relativity, but rather under a subgroup of 
it, manifest in the standard ADM splitting. The local symmetry usually puts constraint 
on the system and only the physical modes appear as propagating modes. However, since 
Horava theory has not full diffeomorphism invariance, one can not obtain full constraint 
to restrict the possible modes to the physical modes. Various aspects of this theory have 
been investigated [20]-[48]. 

Motivated by these, in the present work we are interested in investigating interacting dark 
energy model in the framework of Horava-Lifshitz cosmology. 

2 Horava-Lifshitz cosmology with interacting dark 
energy 

In this section we obtain the equation of state for the dark energy when there is an 
interaction between energy density pd and a dark matter p m . The dynamical variables of 
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Horava-Lifshitz gravity are the lapse and shift functions, iV and JVj respectively, and the 
spatial metric Therefore we can write the metric as: 

ds 2 = -N 2 dt 2 + g i3 [dx l + N l d£)[dx 3 + N j dt), (1) 

The scaling transformation of the coordinates reads (z=3): 

t -> IH and x i -> Zx* . (2) 

Decomposing the gravitational action into a kinetic and a potential part as 

S g = J dtd 3 Xv ^N{C K + C v ) (3) 

and under the assumption of detailed balance [18], which apart form reducing the possible 
terms in the Lagrangian it allows for a quantum inheritance principle [16], the full action 
of Horava-Lifshitz gravity is given by 

r f 2 k 2 k 2 ii f ijk 

S g = j dtd'x^N ~ 2 {K l3 K^ - \K 2 ) - + ^—ifcvX- 



8 13 8(1 - 3A) 



2w*y/g 

i^#2 + AR _ 3A 2 



(4) 



where 



K H = ^ - ViiVj- - V.iV,) , (5) 
is the extrinsic curvature and 

C^' = ^V fc (i2? - -RSI) (6) 

the Cotton tensor, and the covariant derivatives are defined with respect to the spatial 
metric g^. e ljk is the totally antisymmetric unit tensor, A is a dimensionless constant 
and A is a negative constant which is related to the cosmological constant in the IR 
limit. Finally, the variables k, w and \i are constants with mass dimensions —1,0 and 1, 
respectively. 

In order to add the matter component (including both dark and baryonic matter) in the 
theory one can follow two equivalent approaches. The first is to introduce a scalar field 
[20, 21] and thus attribute to dark matter a dynamical behavior, with its energy density 
p m and pressure p m defined through the field kinetic and potential energy. In the second 
approach one adds a cosmological stress-energy tensor to the gravity field equations by 
demanding to recover the usual general relativity formulation in the low-energy limit 
[35, 47]. Similar to [47], where the author add a matter sector to the Horava-Lifshitz 
action, we would like to add matter and dark energy sectors with the following properties: 
It must respect foliated diffeomorphism invariance, obey the principle of detailed balance 
and be nontrivial at the z = 3 critical point and Lorentz invariant in the infrared. So 
inserting a scalar field in the construction and imposing the corresponding symmetries 
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consistently, one results to the following action for the matter field which has interaction 
with dark energy field a (see [20, 21] for non- interacting case): 



S m = dtd 3 x^/gN 



3A - 1 2 m2 



im^V 4 0+^30V 6 0-^(0,a) 



where V((f>,a) acts as a total potential term [49]: 

V t (<j>,a)=V(<j>) + B(<j>,a), 



(7) 
(8) 



the interacting potential B((f),a) is a function ot both fields, and are constants. 

In usual Hofava-Lifshitz cosmological models, the scalar field is responsible for dark 
matter. However, in principle one could include additional scalars in the theory. In this 
work we will allow for one more, in which we attribute the dark energy sector [1]. Thus, 
we add a second scalar a, with action 



S d = J dtd 3 x^N 



^Y^lP + hlh2<jV2(J ~ T^VV + ^<rV 6 <7 - V t ((j), a) 



, (9) 



where hi are constants. Now, in order to focus on cosmological frameworks, we impose 
an flat FRW metric, 



with 



N = l, g l3 = a 2 {t) ll3 , iV* = 0. 



7ijdx l dx j = dr 2 + r 2 d£l\ . 



(10) 



(11) 



We assume that the scalar fields are homogenous, i.e = (pit) and a = a[t). By varying 
N and c/y, we obtain the equations of motion: 



H 2 



K 



6(3A - 1) 



3A- 1 



3n 2 fi 2 A 2 
8(3A - 1) 



+ & 2 )+2V t ( ( p,a)- 



(12) 



3A - 1 



■? + a 2 )-2V t ^a)+ 3KVA2 



8(3A- 1) 



(13) 



where we have defined the Hubble parameter as H = ^. Finally, the equations of motion 
for the scalar fields read: 



(14) 

(15) 



Now we can define the energy density and pressure for the scalar fields. Concerning the 
dark matter, the corresponding relations are 



Pm = 



3A- 1 



Pm = 



<p 2 -v(<p). 



(16) 



Concerning the dark energy sector, we have 

3A - 1 



3/tVA 2 
8(3A- 1) 

3/t 2 /i 2 A 2 
8(3A- 1)' 



Pd = 



& 2 + B((j),a) - 



Pd 



3A-1 . 2 



& 2 -B(<t>,a) + 



(17) 



(18) 



The first parts of these expressions, namely a B((j), a)+ and ^j— a — B(cf),a) cor- 
respond to the energy density and pressure of the a-field, p a and p a respectively. The 
last constant term is just the explicit (negative) cosmological constant. Therefore, in ex- 
pressions (17), (18) we have defined the energy density and pressure for the effective dark 
energy, which incorporates the aforementioned contributions. Using the above definitions, 
we can re-write the Friedmann equations (12), (13) in the standard form: 



H 2 



6(3A - 1) 



Pm + Pde 

2 



K 



H + -II 2 = - 

2 4(3A - r 



Pm + Pde 



(19) 
(20) 



Finally, note that using (14), (15) it is straightforward to see that the aforementioned in- 
teracting dark matter and dark energy quantities verify the following continuity equations 

p m + 3H(p m + p m ) = p m + 3Hp m (l + w m ) = ~^ ~^T~~ = Q ( 21 ) 

p d + 3H( Pd + p d ) = dotp d + 3H Pd (l + w d ) = ^^q^ = -Q- (22) 

The interaction is given by the quantity Q = Tp d . This is a decaying of the dark energy 
component into dark matter with the decay rate T. Taking a ratio of two energy densities 
as r = p m / p d , the above equations lead to 



3Hr 



w d - + 



1 +r r 



3H 



If we define 

pff ^ eff i r 

w d =w d + — , w ™ =Wm ~~3H' 
Then, the continuity equations can be written in their standard form 

p A + 3H(l + wf)p d = 0, 



p rn + 3H(1 + w%)p m = 



Define as usual 



n = — 



K 2 p m 



Q d= P± 

p cr 6(3A-l)iJ 2 ' p cr 

Now we can rewrite the first Friedmann equation as 

fi m + fid = 1. 



« 2 Pd 



6(3A- l)H 2 ' 



(23) 

(24) 

(25) 
(26) 

(27) 
(28) 
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Using Eqs. (27,28) we obtain following relation for ratio of energy densities r as 

r = ^ (29) 

Here as in Ref.[51], we choose the following relation for decay rate 

T = 3b 2 (l+r)H (30) 
with the coupling constant b 2 . Using Eq.(29), the above decay rate take following form 

r _ zb 2 H 

Substitute this relation into Eq.(24), one find 



(31) 



w f = w * + fT d > w ™ =Wm ~T^n~ d (32) 



where 



+ ^0 2 ~V(0) 



Using Eqs. (17), and (27) we obtain following expression 



If a 2 < ~ 12b ^ R2 ; then u^ ff < —1, therefore the interacting dark energy model in the 
framework of Horava gravity exhibiting phantom behavior. 

Now we consider the simplified case of the absence of the a— filed. In this case we have 



-K 2 (B(<t >) - fggffi + 6(3A - l)b 2 H 2 
"-(Wo) - mx _ X) 



Considering the case 3 A - 1 > 0, if B(4>) < f^fzfj , then wf < -1. However in flat 
space non-interacting case, where 6 = 0, wf = w d = 0, so dark energy is a cosmological 
constant and equation of state can not cross over —1. This is big difference with the result 
of [1]. Surprisingly in the absence of the a— filed, when there is an interaction between 
scalar matter filed with cosmological constant we obtain wf < — 1, that is we result to an 
effective dark energy presenting phantom behaviour. This behaviour is pure effect of the 
interaction, so we can ignore the a— filed at first and only allow an interaction between 
matter filed and cosmological constant to obtain the phantom-like behaviour. In this 
case the time evolution of equation of state control by the dynamics of filed 0. 

3 Conclusions 

In order to solve cosmological problems and because the lack of our knowledge, for in- 
stance to determine what could be the best candidate for dark energy to explain the 
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accelerated expansion of universe, the cosmologists try to approach to best results as pre- 
cise as they can by considering all the possibilities they have. Studying the interaction 
between the dark energy and ordinary matter will open a possibility of detecting the dark 
energy. It should be pointed out that evidence was recently provided by the Abell Cluster 
A586 in support of the interaction between dark energy and dark matter [52]. However, 
despite the fact that numerous works have been performed till now, there are no strong 
observational bounds on the strength of this interaction [53]. This weakness to set strin- 
gent (observational or theoretical) constraints on the strength of the coupling between 
dark energy and dark matter stems from our unawareness of the nature and origin of 
dark components of the Universe. May be the recent developments in Horawa gravity 
could offer a dark energy candidate with perhaps better quantum gravitational founda- 
tions. So in the present paper we have formulated Hofava-Lifshitz interacting dark energy 
model. We have considered two scalar fields, one responsible for dark matter fluid and 
one contributing to the dark energy sector. Our calculations show that the interacting 
dark energy model in the framework of Hofava gravity exhibiting phantom behavior. 
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